Tuning the Non-local Spin- Spin Interaction between Quantum Dots with a Magnetic Field 
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1/^ ^ We describe a device where the non-local spin-spin interaction between quantum dots (QDs) can be turned on 

. and off with a small magnetic field. The setup consists of two QDs at the edge of two two-dimensional electron 

' gases (2DEGs). The QDs' spins are coupled through a RKKY-like interaction mediated by the electrons in the 

(N ■ 2DEGs. A magnetic field perpendicular to the plane of the 2DEG is used as a tuning parameter When 

, the cyclotron radius is commensurate with the interdot distance, the spin-spin interaction is amplified by a few 

^ ' orders of magnitude. The sign of the interaction is controlled by finely tuning B^. Our setup allows for several 

' dots to be coupled in a linear arrangement and it is not restricted to nearest-neighbors interaction. 

^ ■ PACS numbers: 73.63-b,73.63.Kv,71.70.Gm,72.15.Qm 
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Quantum information processing requires control and oper- 
ation of interacting quantum mechanical objects UJ. One pos- 
sibility is to produce systems with localized spins in atomic 
impurities, molecules or quantum dots (QDs) and manipulate 
the spin-spin interaction by engineering the electronic wave 
functions of the surrounding material |3l ^ . This would 
allow for the non-local control of spins opening new possibil- 
ities for the fast developing field of spintronics |5|. An impor- 
tant step in this direction was reported very recently by Craig 
et al. P6'l, who coupled two QDs through a confined 2DEG (a 
larger QD) and controlled the magnitude of the interaction by 
closing or opening up the QDs. Besides its relevance for spin- 
tronics, this experiment also opens up the possibility to study 
the interplay between two competing many-body effects: the 
Kondo effect and the RKKY-interaction |0|iH. 

In this work, we propose a different device where the mag- 
nitude (and sign) of the spin-spin interaction between two 
QDs can be tuned by a external field. Our setup consists of 
two QDs at the edge of two semi-infinite two dimensional 
electron gases (2DEGs). When each dot is gated to have an 
odd number of electrons, and therefore a total spin i, they 
interact through the polarization of the 2DEG as in the usual 
case described by the RKKY interaction 1 10]. We show that 
this interaction can be controlled by applying a small mag- 
netic field perpendicular to the plane of the 2DEG. The effect 
relies on the existence of edge states. These are responsible 
for the transverse focusing of electrons injected from a point 
contact [11, 12, 13] or QD |14]. The control mechanism is 
based on the possibility of focusing the electrons that interact 
with one QD onto the other by the action of the external field. 
When the cyclotron radius is commensurate with the inter- 
dot distance, the spin-spin interaction is largely amplified and 
may increase a few orders of magnitude. The mechanism can 
be extended to many spins and to new geometries that allow 
for the independent control of different pair interactions. 

Consider the geometry shown in FiglD It consists of a 
2DEG with a small magnetic field B^^ perpendicular to the 
plane containing the carriers. The 2DEG is divided in two half 
planes by a set of gate contacts which are also used to define 
the QDs as schematically shown in the figure. The Hamilto- 
nian of the system reads 
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FIG. 1: Schematic representation of the proposed device. Both QDs 
are setup to have an odd number of electrons — this is controlled 
by the gate voltages Vi and V2. The QDs spins are then coupled 
through a RKKY-like interaction mediated by the electrons in the 
two 2DEGs. This interaction can be tuned by a perpendicular mag- 
netic field Bz, being maximum when the cyclotron orbit matches the 
distance between the two QD. A third (optional) gate V3 can be used 
to interrupt one of the electrons' path and cancel out the interaction. 



H — Hqd + -ff 2DEG + Hj , 

where the first term is the Hamiltonian of the QDs 
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Here a — 1 and 2 refer to the left and right QD respectively, 
djjg. creates an electron with spin a and energy Ea^r in the QD 
labeled by a and Ua is the Coulomb energy defined by the 
capacitances of the system 1 15]. The single particle energies 
Eacr (measured from the Fermi energy, Ep) can be varied with 
a gate voltage Va as shown in Fig. ^ The second term in 
Hamiltonian Q describes the electrons in the two half planes. 
To describe these 2DEGs, we discretize the space and use a 
tight binding model on a square lattice. 
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where cl^^^ creates an electron with spin a at site n — , Uy) 
of the upper (7 = 1) and lower (7 2) half planes. The 
hopping matrix element t„r„ connects nearest neighbors and 
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includes the effect of the diamagnetic couphng through the 
Peierls substitution. To avoid any zone boundary effects we 
take a lattice parameter qq = 5nm much smaller than the char- 
acteristic Fermi wavelength (Ap ^ 50nm) . We use the Lan- 
dau gauge for which i , = ie~*""^'^'^/'^° and t , ,^.—t, 

where (f) = a^Bz is the magnetic flux per plaquete, (j)Q = hc/e 
is the flux quantum and t — ?i^/2m*a§ with m* is the carriers 
effective mass. The third term in the Hamiltonian describes 
the coupling between the QD and the 2DEGs 
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where c^^^ 



creates an electron at the half 



plane 7 in a linear combination of iVo sites next to the QD a. 
We are interested in the range of magnetic fields that produce 
a cyclotron radius Vc of the order of the interdot distance R, 
defined as the average distance between the sites connected to 
dot 1 and dot 2. For these small fields, the Zeeman splitting 
can be neglected restoring the spin rotational symmetry. In 
what follows we take Ea-^ = Eai- The generalization to the 
case of a large Zeeman energy is straightforward. We assume 
that the QDs are gated to be in the strong Coulomb block- 
ade regime, U + Ea — —Ea, so that the charge fluctuations 
can be eliminated by the Schrieffer-Wolf transformation 
The spin dynamics is then described by a Kondo Hamiltonian 
where Hqo + Hj is replaced by jv^ 



where Sa is the spin operator associated to the QD a and 
1 1 
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with Ta the level width and p the local density of states per 
spin at Ep. and Ua can be measured by transport experi- 
ments llSillQIl . For simplicity, we take ta-y = ta and neglect a 
potential scattering term which is not relevant for the present 
work Finally, usual perturbative procedures give an ef- 
fective exchange interaction between the QDs' spins mediated 
by electrons in the 2DEG. The interdot exchange Hamiltonian 
Hj contains the non-local susceptibility which can be written 
in terms of one-particle propagators. With a negligible Zee- 
man splitting, the propagators are spin independent and the 
Hamiltonian reduces to Hj = J Si ■ S2 with 
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dco/Hlm[Gx(l,2)GT(2,l)], (7) 



where Im denotes the imaginary part and Go- {a, a') = 

{{{ciaa+C2aa), {c\a' a+^la' a))) *e Fouricr transform 
of the retarded Green function L20J and /(cj) is the Fermi 
function. In the following we take ksT ^ J — in this regime 
/(ti>) ~ Q{Ef — Lu). To lowest order in Ja, the one-particle 
propagators are calculated with i/2DEG only. Then, we have 



Gcr{a, a') =gicria, a')+g2aia, a') , 
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FIG. 2: (Color online). Exchange interaction J as a function of the 
perpendicular magnetic field Bz for different interdot distances R = 
0.5 /im (a), 1 fim (b), and 1.5 nm (c) and A'^o — 1 (No = 7 is shown 
in (a) [dashed line]). In the last two cases, J was multiplied by 3 and 
6 respectively. The interaction presents large oscillations when the 
cyclotron radius r-c is such that 2rc — -R/n with n an integer number. 
A fine tuning of allows then to select the magnitude and the sign 
of the exchange interaction. 



with g^a{a,a') = {{c.yac,,c^^a'a))'^+iv- Equation (|7} makes 
evident that, in terms of Feynman diagrams, the effective 
interaction — or the non-local susceptibility — is a bubble di- 
agram with a propagator from dot 1 to dot 2 times a propaga- 
tor from dot 2 to dot 1. In the following, we calculate these 
propagators numerically (see Ref. f\3\ for details). 

The field creates edge states that propagate in oppo- 
site directions on opposite sides of the QDs (see Fig. 0. 
This generates a right-left asymmetry on each half plane, ie. 
g^a (1, 2) ^ g^a (2, 1). In general, with 7^ 0, one of these 
propagators is very small. In fact, with a single half plane, 
the backscattering of the edge states is strongly suppressed 
and they do not contribute to the non-local susceptibility. In 
the proposed geometry, however, tunneling between the two 
half planes described by Hamiltonian Q creates a channel for 
backscattering and the product gi^ (1, 2) g2a (2, 1) does con- 
tribute to the effective interdot coupling. For this reason, an 
effective non-vanishing and controllable spin-spin interaction 
requires a configuration with the two half planes. 

Results for the exchange coupling J are shown in Figure 
12] and |3l normalized by Jq — J{R ^ Xf, = 0). We used 
parameters typical of GaAs systems — an electron density of 
n — 1.5 X 10^^ cm~^, that corresponds to a Fermi energy 
Ef — 5meV measured from the bottom of the conduction 
band and a Fermi momentum fcp ~ 0.1 nm^^. Figure 2 shows 
J as a function of Bz for interdot distances R = 0.5, 1 and 
1.5/im. In all cases, the interaction presents large oscillations 
whenever the magnetic field is such that twice the cyclotron 
radius Tc = TickF/eBz is commensurable with R. We re- 
fer to this fields as the focusing fields since in this situation 
electrons that interact with one dot are focused into the other 
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FIG. 3: (Color online). Exchange interaction as a function of the 
interdot distance R for =0 (a), 231mT (b) and 241mT (c). These 
values of are indicated in Fig. [2^ with arrows. The magnitude 
of J increases a few orders of magnitude when R ~ 2rc — 500 nm 
while its sign depends on the precise value of B^. The solid and 
dashed lines correspond to A^'o = 1 and 7, respectively. 

by the action of the field. As R increases the characteristic 
period of the oscillations and their amplitude decreases (the 
nature of these oscillations is discussed below). The exchange 
coupling J as a function of the interdot distance is shown in 
Fig. |3lfor fixed fields. For = Q the dominant contribu- 
tion is J oc cos(2fcF-R) / (kpR)'^, in contrast to the conventional 
2D behavior J oc cos{2kf: R) / (kpR)^ . The power law decay 
i?^"* is due to the structure of the states near the edge of the 
2DEG. Namely, p depends linearly on energy (it is constant in 
bulk). Since in a semiclassical picture the contribution to the 
propagator Eq. (|8} arises only from the classical trajectories 
near the boundary, one could argue that the effective density 
is oce, so this case "mimics" a 4D one (thus the R^'^ de- 
cay). This is confirmed by both a quantum and a semiclassical 
analytical calculation. For B^^O, a large amplification of J 
is observed for distances such that R — n2rc, where rtis an in- 
teger Comparison of Figs. [3^ andlsji shows that at distances 
of the order of 0.5pm, an increase of the field from zero to its 
focusing value increases the coupling by more than two orders 
of magnitude (notice J/Jq ^ 1 for n — 1). Additionally, the 
sign of the interaction is controlled by finely tuning the mag- 
netic field around the focusing value (see arrows in Fig. [2^). 
Note that a larger Nq (tunneling region) enhances the effect. 
This is due to a reduction of diffraction effects, which leads to 
a better definition of the classical cyclotron orbit 1 13]. 

All these results can be put together in a single density plot 
as shown in Fig. |3 The lines are hyperbolas coiTesponding to 
different focusing fields 

Bi''^ = 2n{hckp/eR) = 2n^Q/XpR. (9) 

Along these lines, ie. when field and distance are simulta- 
neously varied to keep the focusing condition, the exchange 
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FIG. 4: (Color online). Density plot of the exchange interaction J 
as a function of the perpendicular magnetic field B^ and the interdot 
distance R. The solid lines correspond to different focusing condi- 
tions with n — 1, 2, 3 (cf. Eq. Note that the magnitude of the 
interaction is large along those lines while its sign changes (dark and 
light areas). Inset: J as a function of R along the first hyperbola. 
The solid line is a fitting to Eq. JlOt . 

coupling oscillates with large amplitude. A simple fitting of 
the numerical results shows that the dominant contribution to 
J along the hyperbolas is given by 

j^ cos{k,nR/2 + ^) ^ (10) 
R^ 

Note that the interaction decays as R~^ as in the usual 2D 
case — this is consistent with the semiclassical picture since 
now the classical trajectories are not restricted to be close to 
the boundary. The argument of the cosine can be written as 
kpRsff where i?eff = 7r_R/2 is the length of the classical trajec- 
tory (see Fig[3. The period of the oscillations is then equal 
to the period at which the Landau levels cross Ep. In fact, we 
have kpRsff =27ri?eff/AF = 2TTEp/huJc- These oscillations of 
the non-local susceptibility are the analog of the de Haas-van 
Alphen oscillations of the magnetization. The complex oscil- 
lating pattern observed in Fig. [^corresponds to a cut in Fig|3 
along a vertical line. The inset in Figl^shows the coupling J 
when the field and the interdot distance are varied along the 
first hyperbola (n — 1). The solid line is a fitting to Eq. ilQ\ . 

So far we have presented results obtained by fixing the 
chemical potential. However, in a 2DEG with a fixed charge 
density, Ep is pinned at the energy of the partially filled Lan- 
dau level and presents periodic jumps when plotted as a func- 
tion of 1/_B^. Figure|5]shows the coupling J versus the exter- 
nal field for constant electron density at the bulk of the 2DEG. 
Now again, the spin-spin interaction presents a large enhance- 
ment at the focusing fields. There are, however, some differ- 
ences with the previous case: as the external field is varied 
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FIG. 5: Exchange coupling as a function of Bz with the same pa- 
rameters as in Fig. |2|but with the chemical potential pinned at the 
partially filled Landau Level. Note that the sign of the interaction 
tends to be preserved around each focusing field. 



and magnitude. It is important to emphasize that the interac- 
tion between QDs is not restricted to nearest neighbors 

In summary, we showed that two QDs at the edge of a 
2DEG interact with an exchange coupling J that can be con- 
trolled with a small magnetic field perpendicular to the 2DEG. 
When the cyclotron radius rc becomes commensurable with 
the interdot distance R there is a large amplification of the 
interdot interaction. This condition, 2nrc ~ R, defines the 
focusing fields. As the external field is varied around this val- 
ues, the enhanced interaction changes sign allowing for a fine 
tuning of a ferromagnetic or an antiferromagnetic coupUng. 
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3-6343, CNRS-PICS Collaboration Program between France 
and Argentina, and Fundacidn Antorchas, Grant 14169/2L 
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around the focusing values, the sign of the interdot interac- 
tions tends to be preserved. Both ferromagnetic and antifer- 
romagnetic couplings are obtained. The dominant sign of J 
at the different focusing field depends on the parameters, in 
particular on the particle density. The jumps of the Fermi en- 
ergy overestimate some charge redistribution at the edges. In- 
cluding the electron-electron interactions in a self-consistent 
approximation would tend to preserve local charge neutrality. 
This may generate an intermediate situation where the effec- 
tive coupling changes sign as Bz sweeps the focusing values. 

It is worth pointing out that in the presence of a magnetic 
field, the exchange interaction between two spins in the bulk 
of a 2DEG also shows some structure: there is a small en- 
hancement of J for R ~ 2rc while it decays exponentially 
for R > 2rc. In the proposed geometry, however, the focus- 
ing effect produces an amplification of J much larger than 
what is obtained in an homogeneous 2DEG. Moreover, with 
present technologies, it is possible to built a device like the 
one schematically shown in Fig. ^ where the contacts used to 
control the QD parameters {Ea and ta) are also used to divide 
the 2DEG in two halves. An interesting advantage of our setup 
is that it allows to change the magnitude of the exchange cou- 
pling between QDs without modifying their coupling to the 
2DEG. Therefore, transport measurements through the cou- 
pled and decoupled system can be easily compared. 

Since the exchange mechanism requires the two half planes, 
interrupting the particle propagation in one of them decouples 
the QDs. This provides an alternative way to act on the effec- 
tive coupling, which can be implemented with a gate voltage 
on a side contact indicated as V3 in Fig. ^ Also, three or more 
QDs could be built along the central gate with the same or dif- 
ferent interdot distances, allowing for a variety of alternatives 
in which, with the help of side contacts and the external mag- 
netic field Bz, the different couplings could be varied in sign 
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